Annular billiard dynamics in a circularly polarized strong laser field 
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We analyze the dynamics of a valence electron of the buckminsterfullerene molecule (Ceo) sub- 
jected to a circularly polarized laser field by modeling it with the motion of a classical particle in 
an annular billiard. We show that the phase space of the billiard model gives rise to three distinct 
trajectories: "Whispering gallery orbits" , which only hit the outer billiard wall, "daisy orbits" which 
hit both billiard walls (while rotating solely clockwise or counterclockwise for all time), and orbits 
which only visit the downfield part of the billiard, as measured relative to the laser term. These 
trajectories, in general, maintain their distinct features, even as intensity is increased from 10^*^ to 
• cm . We attribute this robust separation of phase space to the existence of twist less tori. 

PACS numbers: 32.80.Rm, 05.45.Ac 



I. INTRODUCTION 

The electrical and chemical properties of fullerenes, 
namely "buckyballs" and nanotubes, remain the focus 
of thorough investigations [1]. The Buckminsterfullerene 
molecule, Ceo, is a prototypical nanocluster because of 
its stability and nearly spherical shape. In particular, 
it is an ideal cage in which to trap so-called endohedral 
atomSjTesulting in molecular systems with peculiar prop- 
erties [2|, i.e., enhanced stability with respect to tempera- 
ture. In recent years, there has been a significant interest 
in subjecting Ceo to extreme conditions to probe its elec- 
tronic and structural stability properties. A new class of 
experiments on Ceo driven by strong laser pulses show 
that its ionization and fragmentation properties are very 
sensitive to the laser intensity and polarization [3]. In 
particular, the yields show remarkable changes with the 
ellipticity of the laser field. 

Motivated by these findings, we consider the motion 
of an electron inside the valence shell of Ceo- The goal 
is to understand the electronic dynamics prior to pho- 
toionization. In strong linearly polarized laser fields, the 
ionized electron can return to the remaining ion by rec- 
olliding with the cage 0, HI when the laser field changes 
sign. This collision can lead to additional ionization or 
even fragmentation Q of the molecule. 

In particular, we investigate the classical dynamics of a 
one-electron model in the shell of fullerene Ceo subjected 
to a strong circularly polarized laser field. In a circu- 
larly polarized field, the dynamics is best visualized in a 
frame co-rotating with the laser field where a conserved 
quantity emerges, the Jacobi constant [6]. All results pre- 
sented in this paper are in the rotating frame in which 
the circularly polarized laser becomes a static field with 
definite upfield and downfield directions. We restrict the 
dynamics to a two-dimensional configuration space (the 
plane of polarization) for the valence electron. The va- 
lence electron feels an averaged potential, which, as we 
later show fSec. Ill A]) , is very close to a spherical square 
well potential, where the electron bounces between the 
walls like a particle in an annular billiard. We choose 



a billiard model for its simplicity both analytically and 
numerically and because it serves as a faithful represen- 
tation of the full model potential while not allowing ion- 
ization to occur. 

Annular billiards occur in the literature in at least two 
contexts: Fermi acceleration and the study of quantum 
chaos by comparison with classical and quantum mechan- 
ical computations [f-T^ . Chaotic dynamics arises either 
from pulsating boundaries or from an off-center inner 
wall. In our treatment the two walls are fixed and con- 
centric. The main distinction here from other works on 
annular billiards is that in the rotating frame the electron 
moves along curved paths between successive wall colli- 
sions. The introduction of a Coriolis term into the Hamil- 
tonian upon the transformation to the rotating frame is, 
of course, akin to introducing an effective uniform mag- 
netic field (and another frequency, the Larmor frequency) 
in which an electron moves on a curved path fl5l-[l9| . The 
laser wavelength is taken as 780nm (corresponding to a 
frequency of 0.0584 a.u.) and its intensity is varied from 
zero to 10^^ W • cm~^, which are values consistent with 
what is routinely performed in experiments on Ceo [!]• 

Our principal finding is that trajectories of the elec- 
tron fall into one of three possible types which originate 
from specific parts of phase space. We identify various 
phase space structures which keep these trajectory types 
distinct from one another. Our classification is as fol- 
lows: First there are "whispering gallery orbits" (WGO). 
These trajectories hit only the outer wall and their di- 
rection of travel, either clockwise or counterclockwise, is 
determined by which of the two regions in phase space 
where they originate. The second type, which we call 
"daisies" from their typical shape, hits both walls and 
can rotate either clockwise or counterclockwise. Both 
daisies and WGOs are called positive or negative based 
upon their direction of travel, counterclockwise or clock- 
wise, respectively. The third type is mainly influenced by 
a very simple elliptic periodic orbit bouncing between the 
two walls in the downfield direction. The shape of these 
trajectories on the Poincare section resembles a popular 
snack food, the Pringles curved potato chip [20|, thus 
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their designation as "pr ingle orbits." 

We show that all the trajectories maintain their dis- 
tinct characteristics with changing intensity because of a 
special class of invariant tori which do not fulfill the usual 
twist hypothesis required by the Kolmogorov-Arnold- 
Moser (KAM) theorem. These tori are usually denoted 
as "twistless" tori or shearless curves (of the Poincare 
map) or "meandering" curves when they are associated 
with separatrix reconnection [2114231 . Numerical studies 
show that these tori are very robust against perturba- 
tion, and are natural candidates for the partitioning of 
phase space into regions where well-defined, qualitatively 
distinct trajectories can be found. 

The plan of the article is as follows: In Sec. |TT1 we 
introduce the model and dynamical rules of the billiard 
and its corresponding Hamiltonian and refiection guide- 
lines. In Sec, nil) we give a qualitative analysis of different 
electron trajectories, and then relate them to the organi- 
zation of phase space. Next, we analyze the dynamics by 
using Poincare section and frequency map analysis |2j 
in order to gain deeper insight into the properties and 
organization of phase space when the laser intensity is in- 
creased, and in particular concerning the role of twistless 
tori in the partitioning of phase space into three principal 
regions. 

II. DYNAMICAL RULES 

A. Hamiltonian Model 

A typical electron inside the shell feels the infiuence 
of the electrostatic potential created both by the posi- 
tive charges from the nucleus and the electronic density, 
combined with the infiuence of the laser field. The effec- 
tive single-particle potential is computed using density 
functional theory from a jellium approximation for the 
positive charge background [25|. It contains steep walls 
in the potential around the radius of fullerene (ro = 6.69 
a.u.) (2514281] with a certain thickness [29]. The Hamilto- 
nian expressed in atomic units (a.u.) and in the dipole 
approximation, reads 

7^(x,p,i) = ^ + y(|x|)+x.E(t) (1) 

where x = (x, y) is the position of the electron in the 
polarization plane, = {Px^Py) its canonically conjugate 
momenta and | • | denotes the Euclidean norm. The circu- 
larly polarized laser field is given by E(t) = £^o(ea? sincjt+ 
e^ycoscjt), where Eq is the electric field amplitude, uo is 
the laser frequency, kept fixed at 0.0584 a.u., and e^^ and 
By are unit vectors along the x and y axes, respectively. 
The laser intensity is the time averaged Poynting vector 
of our laser field and is related to Eq by the relationship 
/ = aE^^ where a = 7.044 x 10^^ when laser intensity is 
measured in W-cm~^. Figure [1] shows the potential F(r), 
where r = |x| as given in Ref. [28]. We note that the po- 
tential is very stiff at the boundaries of the shell. This is 




FIG. 1: Potential V, felt by a valence electron in fullerene as 
given in Ref. 28] . The overhead plane corresponds to the ac- 
cessible billiard region (white space) with a sample trajectory 
(blue curve). 

a common feature of various models for Ceo f25l-[27j. This 
property holds for ions Cgg [26] also. An approximate 
potential consisting of a spherical square well potential, 
where the potential is equal to —Vq for r G]ro — S, rp -|- S[ 
and zero elsewhere, has been proposed in Refs. jsol, Isij . 
This model has succeeded in explaining the oscillations 
in the photoionization cross-section of Ceo [30]. We build 
a billiard model along these lines, where the steep walls 
of the potential are replaced with infinite walls and the 
dynamics in the annular region between the two walls is 
given solely by the interaction of the electron with the 
electric field, later referred to as the laser-driven dynam- 
ics: 

7^(x,p) = M!+x.E(t), (2) 

and refiection rules are applied whenever the trajectory 
reaches r = rin or rout, which are 5.14 a.u. and 8.24 a.u., 
respectively in our computations. 

First, we perform a canonical change of variables into a 
rotating frame (with the laser field). The new coordinates 
{x,y,p^,Py) are given by 

where 

^ ^ y cosut — smut J 

In the new set of variables, the Hamiltonian becomes 
time-independent and reads 

p^ p^ 

JC{x,y,Pa:,Py) = Y + Y -^{xpy - yPa:)^EoX, (3) 

where we have dropped the bars for simplicity. The re- 
sulting Hamiltonian has two degrees of freedom and the 



FIG. 2: Section (y — 0) of the zero- velocity surface in the 
accessible part of the billiard for / = 10^^ W • cm~^ and 
uj — 0.0584 a.u. The red arrows show the deformation of 
the zero velocity surface as intensity is increased. The dashed 
horizontal lines are the Jacobi values used in Fig.[3]and Fig. (8] 



FIG. 3: Accessible regions (white) of the billiard for different 
Jacobi values at / = 10^^ W • cm~^ and cj = 0.0584 a.u. 
Jacobi values correspond to the dashed lines in Fig. [21 



value of the Hamiltonian is the Jacobi constant of celes- 
tial mechanics [6]. 

B. Topology of phase space 

The accessible part of phase space changes depending 
on the field frequency, amplitude, and the value of the 
Jacobi constant. A revealing way to visualize the acces- 
sible part in position space is to compute the zero velocity 
surface [6]. By applying Hamilton's equations to Eq. ([3]) 
we arrive at 



truncations of the annulus are possible; in particular, we 
distinguish two types: one which is homotopic to an an- 
nulus, and one which is only a portion of an annulus (see 
Fig. [3|). For /C larger than + E^oHn, the entire 

annulus is accessible to the dynamics. As the laser in- 
tensity is increased, the difference between the right and 
left sides of the well is amplified. In this paper we mainly 
consider Jacobi constants larger than + E^r-^^ 

such that the full annular region of the billiard is ac- 
cessible to the dynamics. For / = 10^^ W • cm~^ and 
00 = 0.0584 a.u., this critical value of JC is approximately 
equal to 0.15. 



y =Py- ujx, 
so that the Jacobi constant becomes 

/C (X, ^, i, ^) = y + |- - ic^2 ^^2 ^ ^ ^^^^ 

Setting X = y = gives the zero- velocity surface 
V (x, y) = - ^uj^ {x^ + + Eox, 

which charts the lower limit of the accessible parts of 
the billiard, as JC is varied. A cross section of the zero- 
velocity surface is shown in Fig. [2]for = 0. Depending 
on the value of /C, three possibilities arise: If JC is smaller 
than — u'^r^^^/ 2 — Eqt out 1 then there is no accessible part 
to the dynamics. If JC is between —uj'^r'^^^/2 — £^o^out 
and —uj'^rf^/2 + £^onn then only a portion of the annu- 
lar region is accessible. In this range of values, several 



C. Dynamical rules 



The dynamics is computed in a piecewise fashion be- 
cause of the walls. It is composed of segments of laser- 
driven dynamics, as given by Hamiltonian (|3|), until the 
particle reaches one of the walls. At this instant, the re- 
flection rule is applied which mimics an elastic scattering 
at the limit of an infinitely stiff potential. 

Concerning the laser-driven dynamics, the equations 
of motions associated with Hamiltonian (|3]) are given by 
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Px (t) 

Py (t) 



Px,ot I COS out 



Py,o - — ]t^yo ] sincjt, 



^0-^1 -^Px,ot ] sincjt, 



Py,o ) sin out + Px^o COS out, 

Py I COS out — Px sin out H ^ ,(5d) 

ou I ' 



(5a) 

(5b) 
(5c) 



where xq, ^o, Pa;,07 and Py^o are the initial conditions (at 
time t = 0). 

At time t = tn, the electron reaches either one of the 
two walls and we apply the reflection condition before 
a next phase of laser-driven dynamics. The reflection 
conditions are more easily expressed in polar coordinates 
as it corresponds to the radial momentum changing sign 
while the other coordinates are unchanged (see Fig. |4|). 
By definition of the billiard, the rebound takes place at 
time t = tR satisfying 



where i denotes the wall index, i.e., i G {in, out}. Be- 
fore the rebound (i.e., at t = t^), we compute the ra- 
dial momentum pr from the value of {x^y^Px^Py) by 
Pr = {xpx + yPy)/'f' Then, the rebound condition is 
given by a change of sign of the radial momentum 

Pr (4) = -Pr {tn) . 

As a consequence, the values of the momenta after the 
rebound, denoted p^ and p^, are given by 



Pt 



Pt 



(x^ - y^)p^ + 2xypy 
-2xyp- + (x^ - y'^)p- 



(6a) 
(6b) 



We note that Hamiltonian (|3]) is left unchanged by the 
reflection rule, that is JC (t^) = JC (^^), which is easily 
seen from the expression of the kinetic energy in polar 
coordinates which is equal to (2rf). 



D. Linearization of the flow 

In this section, we consider the linear effect of the re- 
bound condition on neighboring trajectories. The moti- 
vation for doing so is twofold: First the tangent rebound 
condition can be used to compute the tangent flow of tra- 
jectories to deduce the linear stability of periodic orbits. 




FIG. 4: Schematic representation for computation of the re- 
bound condition. The solid blue and red lines are real tra- 
jectories, differing only by a small perturbation in the initial 
conditions. The dashed red and blue lines are representa- 
tions of the linearized dynamics. The angle is the angle of 
rotation used in Hamiltonian (|7|. 



Second, as the billiard model corresponds to the limit 
of infinitely stiff potential for a Hamiltonian system, the 
billiard model should preserve the Hamiltonian structure. 
We have already checked that the rebound condition pre- 
serves the Hamiltonian, thus the last prescription is that 
the symplectic two-form is preserved, or equivalently the 
tangent rebound matrix is symplectic. 

The tangent rebound matrix can be seen as an exten- 
sion of the tangent flow [s^ to include the impact of the 
rebound on neighboring trajectories (to first order). For 
that, we consider a first trajectory with initial conditions 
^o^yo^Px,o^Py,o at t = in the neighborhood of one of 
the two walls, i.e., Xq + ?/q ^ r|, where i G {in, out}, and 
such that the rebound time tR <^ 1. Then, we look at 
the impact of a perturbation of these initial conditions to 
xo + dxo, yo + dyo,Px,o + dpx,o,Py,o + dpy^o immediately 
after the rebound. Because the rebound time is not the 
same for the two trajectories, we have to consider a small 
laser-driven propagation before and after the rebound to 
deduce the tangent rebound properties. The situation is 
schematically depicted in Fig. HI 

The computation of the tangent rebound matrix is 
more easily seen in a rotated Cartesian set of coordinates 
x^y^Px^Py, for some angle to be specified later. In the 
rotated frame, the corresponding Hamiltonian reads 



JC{x,y,Px,Py) 



Px , Py / 



ypx) 

-\-Eox cos 6 + Eoy sin 6>, 



(7) 



where we have dropped the tildes for simplicity. Since 
we are only interested in the linear properties of the re- 
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bound condition, throughout this section we consider ah 
the equations hnearized to the first order. 

As explained in Sec. Ill C\ the dynamics is computed 
piecewise with a first stage of laser-driven propagation 
until the electron reaches the wall, then the rebound 
condition and finally a new stage of laser-driven prop- 
agation. In the rotated frame and before the rebound, 
i.e., t < tR^ the laser-driven dynamics given by Hamilto- 
nian ([7]) yields 



x{t) 

Px it) 

Py (t) 



yo + {Py,0 - ^Xo) t, 

Px,0 + {^Py,0 - Eo COS 0) t, 

Py,o - {^Px,o + Eo sin 0) t, 



(8b) 
(8c) 
(8d) 



where we have neglected 0{t^). At this stage, we select 
the rotation angle such that the perpendicular direction 
to the wall at the rebound is aligned with the x-direction 
for the unperturbed trajectory. We keep this fixed frame 
for the perturbed trajectory. An alternative way (which 
provides the same solution) is to consider a perturbed ro- 
tated frame (obtained by a rotation by an angle 0-\-dO) to 
impose the same constraint (the direction perpendicular 
to the wall is the x-axis at the rebound) on the perturbed 
trajectory. As a consequence of the chosen angle 6>, the 
rebound condition solely depends on the x-direction, such 
that X {tR) = Ti, or equivalently using Eq. (|8a|) 



tR^ 



Xo 



Px,0 + ^Vo 



(9) 



Besides, because of the circular shape of the billiard and 
the angle ^, at the rebound the ^/-component vanishes 
{y {tn) = 0). However, for the perturbed trajectory, this 
condition does not apply. Since we consider the per- 
turbed trajectory comparatively to the original one in 
the same rotated frame, it is easier for the purpose of the 
calculation to keep formally the ^/-components, knowing 
that it is actually equal to zero. 

The next step for the trajectory dynamics is the re- 
bound condition, at time t = Ir. Because of the orienta- 
tion of the frame where the x-direction is aligned with the 
radial one, the rebound condition ([6]) becomes x (t^) = 
—X (t^) and y (t^) = y (t^) which implies to the mo- 
menta 



P^ (4) = -P^ {'^r) - (^r) ' 



(10) 



and Py (4) 
changed. 



Py (^i?)' while the positions are left un- 



Finally, after the rebound, the trajectory experiences 
a new phase of laser driven dynamics. Combining the re- 
bound condition ([TO]) with a linearized propagation sim- 
ilar to Eq. ([8j), one can write the dynamics after the re- 
bound as a function of the initial conditions (before the 



rebound), such that 

x{t) ^ x{tR)^{p^{t^)^ujy{tR)){t-tR), 

^ Xo - {px^o + ^^o) {t - 2tR) , (11a) 
y{t) ^ yo^{Py,o-^xo)t, (lib) 
Px {t) ^ - {px,o + 2cj?/o) - 2uj {py^o - OJXq) tR 

+ {ujpy^o - EocosO) it - 2tR) , (11c) 
Py (t) ~ Py,o + 2uj {px,o + ujyo) {t - tR) 

- {ojpx,o Eq sin 0)t, (lid) 

where the equations are linearized at the first order in 
time. With the explicit formula for the dynamics after 
the rebound, as given by Eq. (pT]) . it is straightforward 
to compute the impact of the perturbation on the initial 
conditions by replacing xq, yo^Px,o^Py,o with xo-\-dxo^ yo-\- 
dyo^Px,o + dPx,o^Py,o + dPy,o respectively. Because of the 
change of initial conditions, the rebound time is modified 
to tR + dtR as well. Using Eq. (j9j) for the perturbed 
trajectory we end up with 



dtR 



dxo 



Xo 



■ {dpx,o + (^dyo) 



Px,o + ujyo (p^^o + ^^o) 

(12) 

Finally, combining Eq. (pTj) with Eq. (p!2]) it is possible 
to compute the perturbed dynamics after the rebound. 
We define the deviations dx^dy^ dpx^ dpy of the perturbed 
trajectory after the rebound. For instance, considering 
the x-coordinate, we obtain 



dx 



-dxo 



2 {u - Xo) 



Px,0 + ^^0 

{dpx,o + ojdyo) {t 



{dpx,o + ojdyo) 

2tR). 



(13) 



Since we are interested in the dynamics in the vicinity of 
the rebound, we consider the limits xo Vi and t ^ t^ 
(such that t ^ 0). As a consequence, we end up with 
dx = —dxo from Eq. (p!3|) . A similar procedure can be 
applied to the other components and summarized in the 
Hnear equation 




where Jr is the tangent rebound matrix given by 



Jr 



/ -1 

10 

{Acjpy - 2Eo cos (9 - 2cjVi) p'^ -2uj -1 
V 2uj 1 

(14) 

using the conditions x = , = and where Px is taken 
right before the rebound on the wall. We note that the 
tangent rebound matrix is a symplectic matrix, which 
proves that the rebound condition (j6|) preserves the sym- 
plectic two- form, i.e., 

dxo A dpx,o + dyo A dpy^o = dx A dpx dy A dpy. 
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The tangent flow is used to characterize the hnear sta- 
bihty of invariant structures hke periodic orbits. As for 
trajectories, their integration is carried out piecewise: the 
integration is composed of intervals of laser-driven prop- 
agation and rebound conditions. Between the rebounds, 
we integrate the tangent flow given by [32| 

dtJ = \/F J, 

where VF is the matrix of variations of the flow associ- 
ated with Hamiltonian ([7]). Then, right after a rebound 
on one wall, the Jacobian matrix is equal to the product 
of the previous Jacobian matrix right before the rebound, 
denoted J~, with the rebound matrix (fT4|) : Jacobian ma- 
trix after the rebound reads 



J' 



III. ANALYSIS OF THE DYNAMICS 

A sampling of typical trajectories of the annular bil- 
liard for / = 10^^ W • cm~^ is shown in Fig. [5l These 
examples illustrate qualitatively the different types of ob- 
served trajectories already discussed in Sec. HI In the top 
row we see that the trajectories only hit the outer wall 
and never the inner wall and are therefore "whispering 
gallery orbits", or WGOs. In the middle row we show 
"daisy orbits" . They are qualitatively the same, hitting 
both walls successively and accessing the entire angu- 
lar distribution of the billiard. Both WGO and "daisy 
orbits" keep a constant rotational direction, either clock- 
wise or counterclockwise, which we denote negative or 
positive, respectively. In the bottom row, left panel we 
see a "pringle orbit" , hitting both walls in turn, however, 
limited only to the downfleld region of the billiard. The 
simple two rebound trajectories located at both extremes 
of the upheld and downflled region of the billiard are pe- 
riodic orbits. The leftmost curve is an elliptic periodic 
orbit (stable) while the rightmost curve is a hyperbolic 
periodic orbit (unstable). In the bottom right panel is a 
trajectory which hits neither wall (see Sec. IIII A ip . In 
the following sections we connect these trajectories to 
phase space structures and their stability. 



A. Poincare Sections 

Since the dynamical system has two degrees of free- 
dom, a convenient way to visualize the dynamical orga- 
nization of phase space is by Poincare sections. Here we 
consider a Poincare section with equation = in the 
direction > 0. The rebound condition ([6]) imposes 
a discontinuity in the radial momentum at the re- 
bound, such that formally the condition = is never 
reached during a rebound. However, we see the billiard 
as the limit of an inflnitely stiff potential, and a smoother 
dynamics corresponding to Hamiltonian ([1]) would reach 




FIG. 5: Trajectories of the annular billiard for / = 10^^ W • 
cm~^ and cj = 0.0584 a.u. The trajectory type is shown in 
the top left corner of the panel and IC = 0.35 (see Fig. [2]) for 
all panels with the exception of the bottom right. The arrow 
shows the direction of travel of the trajectory. In the bottom 
left panel the simple two rebound orbits are periodic orbits. 
The downfield orbit (orange) is elliptic and the upfield orbit 
(red) is hyperbolic. The bottom right panel is a trajectory 
which hits neither wall for JC = 2Up ^ 0.002 (see Sec. lIIIAip . 



Pr = before changing sign. As a consequence, we con- 
sider the rebounds on the walls as potential candidates 
for the Poincare section. The electron either rebounds 
on the outer wall, meaning that the radial momentum 
changes from positive to negative values or the electron 
rebounds on the inner wall, so that the radial momentum 
changes from negative to positive values. In order to com- 
ply with the transverse condition p^ > we include only 
the collisions with the inner wall. Furthermore, because 
of the rebound condition (|6]) we note that po^ x = r cos6>, 
and y = rs'mO are continuous under a rebound so that it 
is equivalent to record their values either directly before 
or after the rebound. The Poincare section can be rep- 
resented in several ways. A three-dimensional plot can 
be used where we plot (x^y^po) or we can make a pro- 
jection onto the plane (O^pe) using the constant Jacobi 
constraint. Regardless of which method is chosen there 
are two types of points on the section: The flrst kind are 
points on the inner wall for which r = rin and Pr cho- 
sen so as to satisfy the condition on the Jacobi constant 
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FIG. 6: A Poincare section for / = 10^"^ W-cm-^ cj = 0.0584 
a.u., and AC = 0.35. The blue trajectories are the namesake 
of "pringle orbits". The corresponding two dimensional pro- 
jection is displayed in Fig. [71 bottom right panel. 



/C. The second kind are points for which = and r 
is chosen so as to satisfy JC. In this paper, for the sake 
of simplicity of our figures, we make use of the projec- 
tion onto the two dimensional plane (O^pe)^ however, we 
show the three dimensional counterpart in Fig. [6] where 
we have colored "pringle" trajectories in blue to illustrate 
their namesake shape. 



1. Trajectories which intersect neither wall 

It is natural to ask whether the choice of Poincare sec- 
tion is a good one, i.e., do all trajectories intersect the 
Poincare section? A small subset of trajectories which 
are noteworthy for both their peculiarity and the dynam- 
ics they showcase, do not intersect the Poincare section 
Pr = 0. In that spirit we analyze the trajectories which 
do not collide with either wall. For such trajectories, 
x'^ -\-y'^ needs to remain between rf^ and r^^^ at all times. 
The dynamics of these trajectories is governed by Hamil- 
tonian (j3]). Using the translation x = x — Eq/u'^ , px = Px^ 
y = y and py = Py — Eq/uj^ the Hamiltonian is mapped 
to 



JC 



Px 



- P'^ 

- uj{xpy - ypx) H 



El 



where we have dropped the tildes for simplicity. The 
dynamical features no longer depend on the value of Eq. 
The dynamics is better seen in polar coordinates, where 
the Hamiltonian becomes 



Pr Pe 



El 



Since po is a conserved quantity, the dynamics is that of a 
particle evolving in a potential equal to p^/(2r^) and the 



particle will collide with a wall unless po = {). In the case 
where, = 0, is constant, and it has to vanish so that 
no collision with the walls takes place. Therefore the only 
trajectories which do not hit a wall are circular orbits 
(since r = = 0). In the original coordinates, these 
circular periodic orbits are centered around (xo,yo) = 
(^o/^^^O) and they have a specific Jacobi constant of 
2Up where Up = Eq/{4:Uj'^) is the ponderomotive energy. 
For this Jacobi constant, there exist a priori an infinite 
number of such orbits since the radius is not fixed. The 
only constraint on the radius is that the circular orbit 
has to fit inside the annulus. 

Based on the laser parameters Eq and cj, the exis- 
tence and characterization of such orbits can be divided 
into several categories. For realistic fullerene parame- 
ters, (rout — Hn) /2 < Tin, which is considered here, the 
analysis can be grouped into four categories. Different 
parameters may lead to a different decomposition that 
can nevertheless be identified in a similar fashion (for 
instance the third point below may dissappear). 

1. If Eq/uj'^ is larger than rout, then such trajectories 
do not exist because the center is outside the 
billiard. 



2. If Eo/u^ is in between Tin and Tout? then the 
circular orbits are on the right hand side of the 
inner wall of the annular billiard. 



3. If Eq/uj'^ is in between (rout — Hn) /2 and rin, then 
such trajectories do not exist because none of the 
orbits can be fit inside the allowed region. 



4. If Eq/uj'^ is smaller than (rout — Hn) /2, then the 
circular orbits surround the inner wall with a slight 
shift in the right direction. 



The bottom panel of Fig. [5] shows a sample trajec- 
tory in the fouth category. In this example, the value 
Eo/u'^ ^ 1.1, which is less than rin and hence the tra- 
jectory surrounds the inner wall, but is shifted slightly 
to the right. Except in rare cases (where Eq/uj'^ is equal 
to rin or rout) 7 if such orbits exist, then they exist as a 
continuous family. A linear stability analysis shows that 
these orbits are parabolic. In addition, given that 6 = —u 
(since po =0), the particle turns clockwise. 

Of course, these orbits remain exceptional, in the sense 
that they only exist at some particular value of the Jacobi 
constant JC = 2Up. For uj = 0.0584 a.u., all the orbits hit 
one wall at least for intensities larger than / = 5.56 x 
10^^ W-cm~^ or intensities between 2.16 x 10^^ W-cm~^ 
and 1.96 x 10^^ W-cm~^. The circular orbits confined in 
the right hand side of the annulus only exist for intensities 
between 2.16 x 10^^ W •cm ^ and 5.56 X 10^^ W • cm ^. 
Finally, for intensities lower than 1.96 x 10^^ W • cm~^ 
an infinite family of circular orbits surrounding the inner 
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wall exist. Apart from the examples illustrated in this 
section all trajectories intersect, an infinite number of 
times, the Poincare section. Therefore, we can safely keep 
our definition of the Poincare section without missing 
important dynamics. 



2. Varying the intensity 

For £^0 = 0, Hamiltonian system ([3]) presents a contin- 
uous symmetry by rotation with two degrees of freedom, 
so it is integrable. We show the corresponding Poincare 
section in Fig. [7] (top left panel). As expected, phase 
space is foliated by invariant tori. 

When Eq > the system is no longer integrable and 
some invariant tori are expected to be broken. According 
to KAM theory, a large portion of invariant tori persist 
for £^0 small. In Fig. [3 we show the evolution of phase 
space as the laser intensity is varied. With increasing 
intensity we note the development of a resonance near 
po ~ 2.5. This resonance corresponds to the aforemen- 
tioned two-rebound elliptic periodic orbit shown by the 
left most (downfield) orange curve in the bottom left 
panel of Fig. O This very robust elliptic periodic or- 
bit (situated on the left hand side of the annulus) is 
extremely important in shaping the overall structure of 
phase space as intensity increases. Trajectories originat- 
ing in this region cannot access the entire billiard in a way 
analogous to the librational motion in a pendulum. It is 
this librational motion which yields the "pringle orbits" 
already discussed at the beginning of this section. The 
behavior in the vicinity of this main resonance can be 
described roughly in the following way: The resonance is 
approximately located at po = cjrQ, and this can be seen 
from the dynamical equation for ^, i.e., = —u -\-pe/r'^. 
For a given trajectory, if all the values of pe are larger 
than c^r^y^, then the trajectory turns counterclockwise. 
If the values of po are smaller than ujrf^ then it turns 
clockwise. In between, it oscillates between the two ten- 
dencies. 

The Poincare sections show that the phase space is 
highly regular over several decades of laser intensity. 
Chaotic regions of phase space develop near the hyper- 
bolic periodic orbit which is a rebound between the two 
walls located in the upheld region, or the right-hand side 
of the annulus, and is shown in the rightmost curve in 
the bottom left panel of Fig. O Overall the structure of 
phase space looks very similar to that of a forced pen- 
dulum. In particular, the width of the main resonance 
zone grows like \^Eo or equivalently like /^/^. However, 
we will see that there is a number of discrepancies for 
which twist less tori are the most significant ones. We 
readily observe that the lower part of phase space (nega- 
tive angular momentum) is more chaotic than the upper 
one (positive angular momentum). 
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FIG. 7: Poincare sections for various values of the laser in- 
tensity. Starting in the top left panel we begin with 7 = 
(integrable case) and moving left to right and up to down 
the intensity is increased. In all panels uj — 0.0584 a.u. and 
K = 0.35. 



3. Varying the Jacobi value 

With the intensity fixed at / = 10^^ W • cm~^ the Ja- 
cobi value can also be varied, keeping in mind that this 
variation affects the accessible regions of the billiard (see 
Sec. Ill B|) . We show the corresponding Poincare sections 
in Fig. [HI Starting with the top left panel, the Poincare 
section is contained inside the interval G [2.45, 3.83] 
which agrees with the corresponding panel in Fig. [3l 
note also that the dynamics is highly regular. Besides, 
since the inner wall is not accessible for this Jacobi value, 
points on the section result directly from the condition 
that Pr = 0: A typical trajectory hits the outer wall but 
never reaches the inner one. For the top right panel, 
the inner wall is now accessible (see top right panel of 
Fig. O and the dynamics shows a mixed chaotic and reg- 
ular behavior. In the regular region of the Poincare sec- 
tion the trajectories are regular "pringle orbits" , whereby 
the turning points in are due to the dynamics (and not 
to the geometry of the bilhard). Recall that the regular 
region does not span the entire region of accessible 6 val- 
ues. However, the trajectories originating in the chaotic 
region experience turning points in because they hit the 
artificial walls imposed by the choice of Jacobi value. In 
the bottom panel the entire inner wall is accessible and 
only a small portion of the outer wall is inaccessible. The 
dynamics is still mixed, composed of a regular region with 
"pringle orbits", and a highly chaotic region where the 
trajectories hit the virtual walls imposed by the geome- 
try of the configuration space. WGOs are not possible 
for these values of Jacobi constant. In fact, WGOs ap- 
pear only when the entire annulus is accessible, e.g. in 
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FIG. 8: Poincare sections for various values of AC with / = 
• cm and uj = 0.0584 a.u. These Jacobi values are 
the same used in Fig. [2] and Fig. [3] The top right panel 
includes two trajectories. The left trajectory is taken from 
the "pringle" region and the right trajectory is taken from 
the chaotic region. 



the bottom right panel of Fig. [71 



B. Partitioning of phase space 

Because of the dimensionality of the billiard (two de- 
grees of freedom) invariant KAM tori constitute barriers 
of transport which confine the (chaotic) dynamics to dis- 
tinct regions of phase space. However, they are not ro- 
bust enough to partition phase space at sufficiently high 
value of the intensity. Here, our analysis reveals the ex- 
istence of much more robust invariant objects, namely 
twistless tori which are particularly relevant for the orga- 
nization of phase space [21, 22], since they partition phase 
space at relatively high intensities into regions where the 
different types of trajectories occur. We next introduce 
a diagnostic tool for finding these tori. 



1. Frequency Analysis 

Frequency analysis [2^ is a practical tool in Hamil- 
tonian systems for analysis of the dynamics. For inte- 
grable systems written in action-angle variables (A, (/:?), 
the method consists of plotting the frequency cj(A) = 
dHo/dA as a function of A, which is expected to be 
smooth for (sufficiently smooth) integrable systems. For 
nearly integrable systems, the frequency is computed by 
a windowed Fourier transform of a chosen observable. It 
is computed for an ensemble of trajectories, and plot- 
ted, for instance, as a function of the initial value of 



the action. From this analysis it is possible to identify 
elliptic and hyperbolic islands, regular regions filled by 
KAM tori, and chaotic regions by their respective unique 
signatures. The elliptic islands are expressed as con- 
stant frequency plateaus, the hyperbolic orbits by cusps 
in the frequency, regular regions as apparently continu- 
ous curves, and chaotic regions as non-smooth sections. 
Frequency analysis can also identify regions where the 
twist condition is not satisfied, i.e., when is no longer 
a monotonous function of the action. In this case, the 
standard twist condition for the standard KAM theorem 
is not satisfied, and it gives rise to a new taxonomy of 
dynamical mechanisms, like separatrix reconnect ion and 
twistless tori [2i|,[22|. 

In Figure [9] we plot the frequency as a function of the 
initial momentum po. The ensemble of trajectories are 
the series of points on the Poincare section with initial 
conditions = tt and various po while imposing the Ja- 
cobi constraint to compute the other variables. The mo- 
tivation is that for the integrable case pe is a conserved 
quantity. The frequency analysis of the integrable case 
(Fig.[9l upper panel) shows that the frequency is a contin- 
uously varying function of pe, as is expected. The most 
interesting feature is that the frequency map does not 
change monotonically with pQ . This implies the existence 
of twistless tori located at the extrema of the frequency 
map, dit Po ^ —3.03 and pe ^ 6.11. In the non-integrable 
case (Fig. [9l lower panel), we see the appearence of a 
plateau which corresponds to the "pringle orbits" . In ad- 
dition we notice some chaotic features in the insets, even 
if the overall behavior seems to be quite regular at this 
intensity. Furthermore, the frequency map still exhibits 
two extrema, again which correspond to twistless tori. A 
closer look at the frequency map around po = —3.2 and 
Pe = 6.5 (insets) reveals a rich dynamics with a succes- 
sion of elliptic and hyperbolic orbits, even if the overall 
behavior seems to be quite regular. 



2. Twistless tori 

The aforementioned twistless tori are visualized by 
high resolution Poincare sections. Figure [10] shows 
Poincare sections with initial conditions nearby po ~ 6.5 
(upper panel) and po ~ —3.2 (lower panel), correspond- 
ing to the local maxima in Fig. [9l Both Poincare sections 
are similar in that they show well-developed chaotic re- 
gions sandwiched between two regular regions. The reg- 
ular regions correspond to WGOs, shown in blue, and 
"daisy orbits" , shown in red. Likewise, the chaotic region 
in both panels exhibits a stratification whereby a trajec- 
tory originating in one of these regions remains there and 
cannot pass to another chaotic region. The stratification 
is not due to KAM tori since KAM tori come in fam- 
ilies. Instead it is caused by the existence of twistless 
tori, which, having dimension two, can partition phase 
space. In particular in the lower panel of Fig. [TOl we rec- 
ognize one of the signatures of twistless tori which is the 
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where each chaotic region exists above (below) the black 
line gives the accessible region of the inner wall to the 
electron. Figure [TT] illustrates how the accessible region 
of the inner wall changes as one moves from WGOs to 
"daisies" in the lower panel of Fig. [TOl Beginning with 
the WGOs, shown in blue, which have no points above 
the black line (no accessible region on the inner wall) we 
increase pQ to the light blue region there is a small range 
of values for which the chaotic region is above the black 
line. The values are centered about ^ = and they 
correspond to the area along the inner wall which is ac- 
cessible to the electron. Moving again upwards in po we 
pass by several more chaotic regions, each having a wider 
range of 6 above the black line which allows for a wider 
range of the inner wall to be visited by the electron. The 
inner wall becomes more accessible until finally reaching 
the "daisy" region, shown in red. At this point, all inter- 
sections of the Poincare section are above the black line 
covering ^ G [0, 27r] and therefore the entire inner wall 
is accessible. The four twistless tori which exist in this 
chaotic region are responsible for the discrete transition 
from WGOs to "daisy orbits" . Likewise, a simiar feature 
can be observed for the upper panel of Fig. [TOl nearby 
pe = 6.5, however there are fewer chaotic regions and 
hence fewer twistless tori. 



FIG. 9: Frequency analysis in the integrable case (/ = 0, 
upper panel) and non- integrable case (/ = 10^^ W • cm~^, 
lower panel). The insets display the frequency map around 
the extrema. In both panels JC = 0.35 and u = 0.0584 a.u. 



meandering behavior (see for instance the interface be- 
tween the orange and the dark red chaotic regions). The 
black line, superimposed over both panels, separates the 
two different possible ways of intersecting the Poincare 
section. Points on the section above (resp. below) the 
black line for the upper (lower) panel are standard in- 
tersections of the fiow with the Poincare section in the 
sense that Pr changes sign smoothly before and after the 
intersection. Points below (resp. above) the black curve 
for the upper (lower) panel are collisions with the inner 
wall where the sign of pr changes due to the rebound 
condition (see Sec. IIII Ap . As expected, all the WGO 
trajectories are below the black curve in the lower panel 
(since none of their points intersect the inner wall). The 
entire WGO region is regular. We also notice that all 
the "daisy orbits" are above the black curve, and this 
region is also mostly regular. Between these two regions 
is a strongly chaotic region with very few elliptic islands. 
Each of the chaotic regions (in both panels) has a por- 
tion above and below the black line. The range of 



Conclusions and outlook 



The motivation of our work stems from recent pho- 
toionization experiments in strong field of atoms and 
molecules in both circular and linear polarized light 
where a significant variation of the yields with polar- 
ization was observed [3]. We propose a rather simple 
dynamical model for the motion of a valence electron in- 
side the valence shell of fullerene Ceo, namely an annular 
billiard. We have investigated the dynamics when this 
electron is subjected to a circularly polarized laser field. 
We have shown that it exhibits three distinct types of 
trajectories: "whispering gallery" , "daisy" and "pringle" 
orbits. These trajectories are found in distinct, identifi- 
able regions of phase space for a wide range of laser inten- 
sity and Jacobi values. They are kept characteristically 
segregated from each other by the existence of twistless 
tori which partition phase space. The twistless tori are 
identified through a frequency analysis and are confirmed 
by generation of high resolution Poincare sections. These 
twistless tori, blessed with high stability, exist in chaotic 
regions where KAM tori have been broken by the strong 
laser field. Because of the barriers they create, twistless 
tori, allow for a transition scenario from WGO orbits to 
"daisy orbits" in both rotational directions, positive and 
negative. 
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